Elastic moduli of vortex lattices within nonlocal London model 
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Vortex lattice (VL) elastic response is analyzed within nonlocal London model. The squash 
modulus turns zero at the field Ha where VL undergoes square-to-rhombus structural transition. 
In the field region H > Ha where the square VL is stable, the rotation modulus turns zero at a field 
H r indicating instability of the square VL with respect to rotations. The shear modulus depends on 
the shear direction; the dependence is strong in the vicinity of Hu where the square VL is soft with 
respect to the shear along [110]. The H dependences of elastic moduli are evaluated for LuNi2B2C. 

PACS numbers: 74.60.-w, 74.60.Ge, 74.70.Dd 
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The theory of elasticity for vortex lattices (VL) in type- 
II superconductors is commonly constructed in a man- 
ner similar to the standard elasticity |Q . One introduces 
the local VL displacement u from equilibrium and ex- 
pands the elastic energy E in a series of derivatives Ui : k = 
dui/dxk- The similarity, however, ends already at this 
stage. Elastic energy of an element of a solid does not de- 
pend on rigid rotations of this element, so that E depends 
only on symmetric combinations uik = (ui_k + Uk,i)/2 
which form the strain tensor, whereas the antisymmetric 
combinations tOik = {ui.k — Uk,i)/2 representing rotations 
are irrelevant. This property can be tracked back to the 
basic isotropy of the space in which solids are situated. 

Vortices, however, are hosted by crystals which at best 
have the cubic symmetry. In other words, the "space 
in which vortices live" is never isotropic. The crys- 
tal anisotropy affects intervortex interactions in a trivial 
way through the tensor of "superconducting masses" 
which determines the anisotropy of the London penetra- 
tion depth A and the coherence length £. For all sym- 
metries except cubic (for which rriik reduces to Sik), rigid 
rotations of VL cause an energy change even for rotations 
about the direction of vortex axes (that of the magnetic 
induction B) @j. 

Moreover, the current and field distributions around 
the vortex core are affected by the underlying crystal in 
a more subtle manner than via the A anisotropy. In full, 
this influence can be described only by the microscopic 
theory. In materials with A/£ = K 1 the problem is 
simplified at distances r>( relevant for intervortex in- 
teractions in fields well below the upper critical field H C 2- 
This is done in the frame of the London theory corrected 
for the nonlocality of the relation between the current and 
the vector potential. The kernel of this nonlocal relation 
provides the formal bridge between the Fermi system of 
electrons in a given crystal and interacting vortices in the 
superconducting condensate ||,|| . In general, this kernel 
depends on the Fermi surface and on the symmetry of 
the order parameter [pi; for simplicity, we consider here 
the isotropic order parameter, so that all superconduct- 



ing anisotropies are determined by the Fermi surface. We 
focus on the simplest part of the VL elastic response and 
consider only straight and parallel vortices thus restrict- 
ing ourselves to the two-dimensional (2D) problems. In 
fact, we calculate only "shear" moduli (meaning defor- 
mations preserving the vortex density, divu = u iti = 0). 
These moduli are practically non-dispersive, that allows 
us to consider uniform deformations (k = 0) Q . 

We begin with the general discussion of the deformed 
VL in cubic or tetragonal crystal in fields along the 
crystal direction c, which allows us to correctly count 
independent clastic constants. Then we evaluate the 
moduli using the London model with nonlocal correc- 
tions for LUNI2B2C. As expected, the squash modulus, 
corresponding to the transformation of the rhombic unit 
cell to the square, turns zero at the transition field Hu 
0. We find also novel features of the elastic response. 
The shear moduli depend strongly on the displacement 
direction, the shear "polarization" : the shear polarized 
along [110] is much softer than others near H D . The 
rotation modulus of the square VL turns zero at a field 
H r > H a , indicating instability of certain orientations 
of the square VL. Then we discuss implications of our 
results for the VL physics. 

The general form of the elastic energy in terms 
of strains and rotations is E = {\jkiUijUki + 
r]ijkiUijWki)/2 + QjkiUijUJki ]§]■ This form is not con- 
venient for counting number of independent coefficients 
because one has to consider symmetry restrictions upon 
three 4th rank tensors. A more direct approach is to 
forgo the splitting of u^k into symmetric and antisym- 
metric parts (which no longer simplifies the problem) and 
to write the energy in terms of itj 



E — ■~fikl m Ui ! kUl,m/2 ■ 



(1) 



The tensor jikim is not symmetric relative to i <-> k and 
I <-> m, but preserves the property of the standard elastic- 
ity tensor 'jikjm = lim,ik- The symmetries of this tensor 
and the number of its independent components are de- 
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termined both by the crystal and by the equilibrium VL 
structure. In other words, "fikim has only the symmetries 
common to the crystal and the equilibrium VL. 

The equilibrium VL in cubic or tetragonal crystals in 
fields above H a along [001] has a square unit cell with 
sides along [110] and [110] §-|l|,?]. We choose the co- 
ordinates x,y in the ab plane as [100] and [010]. Then, 
x may enter the indices of 7 only even number of times, 
and the components obtained from each other by replace- 
ments x «-> y are equal: 

E ~ [^/xxxxi^x^x ^y.y) ^^fxxyy^x,xUy,y 

~t~ ^xyxyi^x^y ^y,x) ^Ixyyx^x^y ^y ,x\ ■ (2) 

After excluding compressions (B is constant, divu = 
Ux,x + Uy.y = 0) we are left with three constants: 

E = (Tim — lxxyy) u x ,x + ~^lxyxy(u xy + u y,x) 
"T ^fxyyx^x ,y^y ,x 

= J [7i u l,x + 72 {ul iV + u 2 y x ) + 2j s u Xty u yiX \ . (3) 

How to classify the three constants, to a large extent, 
is a question of semantics. For our particular problem, a 
uniform deformation defined as 

u = n(xe x - ye y ) , (4) 

(/i is a small constant and e x ,e y are unit vectors) is of 
a special interest because it is this deformation which 
transforms the square VL above H a into a rhombic 
one for H < H a - Since this is a second order phase 
transition, the corresponding modulus must vanish at 
H = H a • This deformation was named "squash" Q . 
For the displacement (|4|), u XiX — —u y y — /i, whereas 

^y,x V'x 

. y — 0. Then, the part of the energy re- 
lated to the squash is E sq = ji[i 2 /2. In other words, 

71 = 2(-y xxxx — Jxxyy) can be called the squash modulus 
and denoted as C sq . 

Consider now a uniform shear polarized along x, 

u = uye x , (5) 

for which only u x , y ^ 0. The energy is E x = j2U x y /2 = 
72^ 2 /2. The energy is the same for the y polarization: 
E y = E x = j2V 2 /2. Hence, 72 is the shear modu- 
lus for either x or y polarizations; below we comple- 
ment the traditional notation with polarization index: 

72 = C66,x = Ce6,y 

Let now the shear be directed along some direction e x i 
so that the displacement is given by u = v y ' e x > , where 

x' = x cos 6 -\- y sin 8 , y' — — xsind + ycos8 (6) 

and 8 is the angle between x ' and x. In the old coordi- 
nates, 



u = v (—x sin 8 + y cos 8) (e x cos 8 + e y sin 8) . (7) 

Then u x>x = —u y . y = — v sin8 cos8, u XiV = vcos 2 8, 
Uy.x = — sin 2 8, and 

v 2 

E{8) = y[72 + (71 - 272 - 273) sin 2 8 cos 2 8] . (8) 

In particular, E(tt/4) = v 2 {~/i + 272 - 2^)j%. Since all 
uf j — ^ 2 /4, the factor 71 + 2(72 — 73) is the modulus 
for the shear polarized along [110]. Thus, in addition to 
C sq = 71 , there are two independent shear moduli in the 
problem: CWx = 72 and C 66 (ir/4) = 71 + 2(72 - 73). 
Consider now a rigid rotation 

u = uj(xe y - ye x ) , (9) 

for which u y>x — —u x ^ y — lu, whereas u x _ x = 0. Then, 
E r = uj 2 (j2 — 73) = C r uj 2 /2. In other words, the ro- 
tational modulus is C r = 2(72 — 73). Thus, the three 
independent moduli can also be chosen as 

C sq = 71, C e6 . x = 72, and C r = 2(72 - 73). (10) 

Note that 

C 66 (7r/4) = C sq + C r . (11) 

The choice of a particular set of elastic moduli is a matter 
of convenience in a problem at hand. 

The above enumeration of the elastic constants per- 
tains only for H > H n (parallel to c) where both the 
crystal and the VL have the square symmetry in the xy 
plane. In fields under H n , the equilibrium VL is rhom- 
bic, and the replacement x «-» y is no longer a symmetry 
operation. We then obtain instead of Eq. (||): 

B ~ i^fxxxx^x.x ^iyyyy^y ,y ^'fxxyy'^x.x'^y^y 

+ lxyxyU x ^ y + JyxyxU y x ) + JxyyxU x ,yUy^ x . (12) 

For an incompressible VL, we obtain four independent 
constants. We can choose them as the squash C sq , two 
shears 6*66.2; 7^ Cee. y , and the rotation C r . 

We now turn to evaluation of the moduli using the 
London equations corrected for nonlocality, which are 
obtained from the BCS theory for superconductors with 
isotropic gap ||: 

— ji = --72 ("ly 1 - A 2 ?iy7 m fcifc m )aj . (13) 

Here, my is the normalized mass tensor (detm^ = 1) 
and A is the average penetration depth (A 3 = A1A2A3, 
Ai = A^/mi, ...) Further, a = A + 4>^V8/2tt and 0o is 
flux quantum. Tensor riiji m is defined as 
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ZTi 2 (v l v.jV l v m ) 
4(?; 2 )A 2 A 2 



E/3- 2 /3i _1 



(14) 

where v is the Fermi velocity and the brackets (. . .) stand 
for averages over the Fermi surface; [3 2 = A 2 + h 2 u) 2 , 
A(T) is the uniform BCS gap, %uj = 7rT(2n + 1) with an 
integer n, and sums run over u> > 0; f}\ — [3 + h/2r with 
r being the scattering time due to nonmagnetic impuri- 
ties. In dirty superconductors 7 ~ (tA/H) 2 — > 0, i.e., 
nonlocal effects vanish. In the clean limit 7 ranges from 
2/3 at T = to w 0.3 at T = T c . 

The flux quantization V x a = h— <fio e z S(r) combined 
with Eq. ([[3]) gives the field of a single vortex along c of 
cubic or tetragonal materials: 



Mk) = 



l + A 2 fc 2 + A 4 K^fc 4 + dfc 2 fc 2 ) ' 



(15) 



where d = 2n 7 



6n 



xxyy- 



The free energy density 



is a sum of pairwise interactions of vortices jljj ] ; in the 
reciprocal space we have 



F 



B 2 

87T0Q 



£%(q)/i,(q) 



(16) 



where q form the reciprocal VL. The factor s(q) is intro- 
duced to deal with the shortcoming of the London model 
which disregards spatial variations of the order parameter 
in vortex cores. At low temperatures of our interest, we 
use s(q) = exp(— q 2 £ 2 /2) [|1|. By and large, this form of 
the cutoff is confirmed in neutron scattering experiments 
from which the form- factor \s(q)h z (q)\ 2 can be extracted 
p6| . We consider the typical case of a platelet sample 
with the c axis normal to the flat face in fields along c. 
Then minimization of F at a given B yields the equilib- 
rium VL, whereas its elastic properties are evaluated by 
considering deviations from the minimum. 

For high-K materials, away of the lower critical field, 
the interaction of vortices extends over several intervor- 
tex spacings. This gives rise to the nonlocal elastic re- 
sponse of VL's ||, for which the elastic moduli, in gen- 
eral, depend on the deformation length scale (on the wave 
vector k). However, even in isotropic case the shear mod- 
ulus is practically k independent. We focus here on the 
nondispersive modes, i.e., we consider the elastic response 
to deformations with k = 0. 

We begin with the equilibrium square VL with basis 
vectors a! = a (e x + e y )/\/2, a 2 = a (e y - e x )/y/2. 
The reciprocal lattice is 



7iV2(m + n)/a , 



TT\2{m — n)/ao , (17) 



with integers m and n. For the 2D uniform deformations, 
the displacement with constant Ui j 's. Since 

the VL cell area is fixed at a given £?, u x>x +u y<y — [jTj. 
The deformed cell is given by a' ai = a a i + Uija a j with 



a = 1,2. It is readily shown that the reciprocal VL is 
also homogeneously deformed q' = q + u* with 



Ua = — u 



(18) 



(with the tensor of derivatives being transposed). 

We now expand the intervortex interaction V(q') = 
s(q')h z (q')/0o i n powers of u*: 

vM-v W+ «££4«.;*k + .... (19) 

The term linear in u* in the sum ( |T6| ) must vanish since 
u* = correspond to equilibrium. We then obtain for 
the elastic energy E — -F{q'} — 

E =T^^{dq-^ q ^h^- (20) 



For the squash deformation (0) this gives 



B 2 ( 2d 2 V 2 d 2 V d 2 V , 

Cs i = ^2^\ + - 2 ^iy7nr^r ) ■ ( 21 ) 



lX dql dql 



dq x dq y 



Similarly, we obtain for the shear (||) and rotation @; 

(22) 



B 2 ^ 2 d 2 v 



>!l 



B 2 ^ ( ,d 2 V 2 d 2 V d 2 V , 

% M +qx H~ qxqy d^dq-y)- (23) 



These three moduli form the complete set of 2D inde- 
pendent elastic constants for the square VL in cubic or 
tetragonal crystals in fields H > H n along [001]. If 
needed, C^(-k/A) can be evaluated using Eq. (pi]). 

In a similar manner one can obtain the four moduli 
for the rhombic VL below H a . The moduli can be cal- 
culated numerically by evaluating the sums in Eqs. ( |2l| 
- ^3|). The results of such a calculation for LUNI2B2C 
(with parameters given in Ref. [El) are shown in Fig.l, 
where the field dependence of the moduli is shown. It is 
worth noting that C r turns zero at some field H r which is 
about twice as big as Hn ■ This implies instability of the 
square VL with diagonals along a and b (which is stable 
for H a < H < H r ) with respect to rotations. For fields 
H > H r we have four equilibrium square VL's of the 
same energy: the two with diagonals rotated relative to 
[100] counter- and clockwise over an angle if, and two for 
rotations relative to [010]. We have found by the direct 
numerical energy minimization that ip = at H = H r 
and increases as \JH — H r for (H — H r ) <C H r . For 
LuNi 2 B 2 C we estimate that ip may reach 7-^8 degrees. 
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FIG. 1. The shear, rotation, and squash moduli versus field. 




FIG. 2. The ratio C 66 ,x /C B6 (tt /4) versus field. 

As Fig.l shows, Ce6,x increases with field. For 
H < Ha, C66,x 7^ C66,y, whereas in the square phase 
between H a and H r they coincide, as they should by sym- 
metry (the moduli for B ~ 200 G are not plotted because 
the VL in domain differs from that for 200 G< B < H n ). 
However, response to the shear other than of x or y polar- 
ization, differs from either Cqq <x or Cqq^. The difference 
is maximum for the shear along [110] in the vicinity of 
the transition field H a , see Fig. 2. 

We thus come to a striking conclusion: weak nonlo- 
cal corrections to the London description suffice to make 
elastic properties of VL's strongly anisotropic even in the 
ab plane where cubic or tetragonal materials are macro- 
scopically isotropic. 

Softening of the squash mode near H a should affect 
vortex fluctuations. This, however, does not lead to ei- 



ther fluctuations divergence or to a peak in the critical 
current Q , the question recently discussed in Ref . (ijj . 

The rotation modulus C r vanishing at H — H r affects 
the equilibrium VL structure above H r . Although some 
indications for the square VL being unstable relative to 
rotations in sufficiently high fields were seen in decora- 
tions of LuNi2B2C |l^], experimental confirmation of this 
instability is still to be established. 

The square VL is soft with respect to shear displace- 
ment along the square sides (along [110] or [110]). This 
suggests that the vortex rows parallel to one of these di- 
rections can easily slide by each other. If this happens 
for, e.g., [110] direction, the VL will possess a long range 
order along [110], but will look disordered if observed 
along [110]. Such situations have been seen in decora- 
tions in fields near H a jig] . The anisotropy of the shear 
moduli can, in principle, be probed experimentally if the 
critical current is arranged to flow only along [100] or 
[110] directions. 

Ames Laboratory is operated for US DOE by the Iowa 
State University under Contract No. W-7405-Eng-82. 
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